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Abstract 

The classical notion of a centralized clock that governs circadian rhythmicity has been challenged with the discovery of 
peripheral oscillators that enable organisms to cope with daily changes in their environment. The present study aimed to 
identify the molecular clock components in Atlantic cod {Gadus morhua) and to investigate their daily gene expression in 
fast skeletal muscle. Atlantic cod clock genes were closely related to their orthologs in teleosts and tetrapods. Synteny was 
conserved to varying degrees in the majority of the 18 clock genes examined. In particular, aryl hydrocarbon receptor nuclear 
translocator-like 2 {arntl2), RAR-related orphan receptor A (rora) and timeless (tim) displayed high degrees of conservation. 
Expression profiling during the early ontogenesis revealed that some transcripts were maternally transferred, namely arntH, 
cryptochrome lb and 2 (crylb and cry2), and period 2a and 2b {per2a and per2b). Most clock genes were ubiquitously 
expressed in various tissues, suggesting the possible existence of multiple peripheral clock systems in Atlantic cod. In 
particular, they were all detected in fast skeletal muscle, with the exception of neuronal PAS (Per-Arnt-Single-minded) 
domain-containing protein (npasi) and rora. Rhythmicity analysis revealed 8 clock genes with daily rhythmic expression, 
namely arntl2, circadian locomotor output cycles kaput {clock), npas2, cry2, cry3 per2a, nuclear receptor subfamily 1, group 
D, member 1 (nridl), and nrld2a. Transcript levels of the myogenic genes myogenic factor 5 (myfS) and muscleblind-like 1 
(mbnll) strongly correlated with clock gene expression. This is the first study to unravel the molecular components of 
peripheral clocks in Atlantic cod. Taken together, our data suggest that the putative clock system in fast skeletal muscle of 
Atlantic cod has regulatory implications on muscle physiology, particularly in the expression of genes related to 
myogenesis. 
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Introduction 

Circadian clocks are ubiquitous time-keeping systems found in a 
wide range of organisms, enabling them to adjust to light (L) - dark 
(D) cycles over a 24 h period [1]. It is not surprising that circadian 
clocks are evolutionarily conserved [2,3] since synchronization of 
behavioral, physiological and metabolic processes to cyclic 
environmental factors represent a significant advantage [4]. The 
core machinery of the circadian clock consists of intracellular 
transcriptional-translational feedback loops regulated by clock 
genes and corresponding proteins [1]. This regulatory network is 
composed of two main arms: the transcriptional activation 
(positive) and the repressor (negative feedback) arms [5]. The 
positive arm of the core clock is formed by two members: circadian 
locomotor output cycles kaput (Clock) and Aryl hydrocarbon 
receptor nuclear translocator-like (Arndl or Bmall). They are 
basic-helix-loop-helix (bHLH) transcription factors that hetero- 



dimerize via a PAS domain in the cytosol [6]. The Clock: Arntll 
heterodimer binds to E box (CACGTG) enhancer elements in the 
promoters of genes encoding period (per) and cryptochrome [ay) in the 
negative arm of the transcriptional-translational loop. In the 
negative feedback loop, Per and Cry proteins block the action of 
Clock: Arntll and thereby inhibit their transcription. However, in 
some organisms either Clock or Npas2 proteins form heterodimers 
with Arntl proteins to activate transcription (Clock: Arntll and 
Npas2:Arntll) [7]. In Drosophila, TIM and PER form heterodimers 
that act as negative regulators of the CLOCK-BMAL1 complex 
[8]. In the stabilizing loop, Clock:Arntll activates the orphan 
nuclear receptors Rev-erb ot/p and RAR-related orphan receptor 
alpha (Rora). These genes are part of the core clock, as they link 
the feedback loops through their proteins that repress or activate 
arntll expression. This autoregulatory mechanism results in a 
cyclic, self-sustained expression of clock genes with an approxi- 
mately 24-h period [1]. 
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In teleosts, clock genes have been studied mainly in model 
species such as zebrafish (Danio rerio), tiger pufferfish (Takifugu 
rubripes), spotted green pufferfish (Tetraodon nigrovirides), Japanese 
medaka [Oryzias latipes) and the three-spined stickleback (Gasterosteus 
aculeatus) [1,23,9]. The molecular components of the clock system 
in a number of aquaculture fish species, such as in European sea 
bass (Dicentrarchus labrax) [10], Senegalese sole (Solea senegalm.su) 
[1 1], rainbow trout (Oncorhynchus mykiss) [12] and Atlantic salmon 
(Salmo salar) [13], have also been studied to a limited extent, he 
Somalian cavefish (Phreatichthys andruzzii) provided significant new 
insights into the evolution and entrainment of circadian clocks. In 
spite of being blind, this cavefish still retains a clock that responds 
to food rather than light [14]. A recent study in medaka has 
demonstrated that the onset of clock rhythmicity can be 
manipulated by accelerating development through elevated 
temperature or by artificially removing the choriont [9]. 

The central clock in teleosts is believed to be located in the 
pineal gland (or retina) [15]. However, clock genes are not only 
confined to these central pacemakers but they are also present in 
various tissues, thus implying the possible existence of peripheral 
circadian oscillators [16]. Circadian oscillators discovered in 
peripheral tissues challenge the hypothesis that circadian rhyth- 
micity is controlled exclusively by the central clock machinery. In 
fact, it has been demonstrated that peripheral clocks have an 
integral role in a specific tissue or cell as they drive the circadian 
expression of specific genes involved in various physiological 
processes [17]. The study of peripheral oscillators has been of 
particular significance in fish because some peripheral tissues are 
photosensitive and the circadian clocks present within individual 
cells can be directly entrained by light [18]. In spite of their 
importance in driving the circadian expression of specific genes 
involved in muscle physiology [6,19,20], peripheral clocks in 
skeletal muscle have received scant attention in fish, with the 
exception of a single study in zebrafish [21]. Adantic cod (Gadus 
morhud) is a particularly interesting species to study circadian 
rhythms in skeletal muscle, since their growth is markedly 
influenced by photoperiod conditions. Continuous illumination 
significandy promotes body growth [22] and modulates the 
expression of myogenic regulatory factors in fast skeletal muscle 
of Atlantic cod [23]. The existence of a peripheral clock in Atlantic 
cod skeletal muscle is plausible and it could well influence 
myogenesis. It has been shown in mice that myoD, a well-known 
skeletal muscle-specific transcription factor which regulates the 
myogenic program, is directly regulated by Clock and Email [19]. 

The goal of this study was to characterize the molecular 
components of the circadian clock system in Atlantic cod fast 
skeletal muscle. Eighteen putative clock genes have been cloned 
and characterized in this study. Their expression was determined 
during a daily cycle (12L:12D) and correlated to the expression of 
muscle-related genes. 

Materials and Methods 

Animal ethics 

All fish handling procedures in this experiment were in 
accordance to the guidelines set by the National Animal Research 
Authority (Forsoksdyrutvalget, Norway) and have been approved 
by the ethics committee of the Faculty of Biosciences and 
Aquaculture of the University of Nordland, Norway. 

Molecular cloning of Atlantic cod clock genes 

Previously characterized clock genes from several species were 
used to identify the clock genes of Atlantic cod by homology 
cloning. To amplify the cDNA fragments, gene specific primers 



were designed either from contigs of the assembled expressed 
sequence tags obtained from the National Center for Biotechnol- 
ogy Information (NCBI, www.ncbi.nlm-nih.gov) database or the 
Ensembl Atlantic cod genome assembly (www.ensembl.org; 
gadMorl). 

Partial gene sequences were amplified using pooled cDNA 
(brain, pituitary and muscle) derived from three adult individuals 
(sex undetermined) following a PGR amplification protocol 
involving 40 cycles of denaturation (30 s at 94°C), annealing 
(30 s at 60°C), and extension (1 min at 72°C). The resulting PGR 
products were cloned and sequenced as detailed elsewhere [24]. 
The identity of the amplified cDNA fragment was confirmed by 
BLAST similarity searches at NCBI. 

Phylogenetic and synteny analyses of Atlantic cod clock 
genes 

Partial sequences of Adantic cod clock genes were used as 
probes for in silico cloning using the Atlantic cod genome assembly 
available at Ensembl, in order to obtain longer sequences for 
phylogenetic analysis. Deduced amino acid sequences were 
aligned with their corresponding orthologs in various species 
using MUSCLE Multiple Sequence Alignment (www.ebi.ac.uk/ 
Tools/msa/muscle). The aligned sequences were further fine- 
tuned by eliminating gaps and highly divergent regions with 
Gblocks 0.91b (molevol.cmima.csic.es). The resulting trimmed 
aligned sequences were used for phylogenetic analysis in PhyML 
3.0 (www.atgc-montpellier.fr/phyml). Likelihood analysis was 
performed using the LG substitution model with four substitution 
rate categories and an estimated y shape parameter. Branch 
support was calculated by aLRT SH-like tests. TreeDyn 
(TreeDyn.org) was used to view the resulting radiation trees. 
Synteny analyses were performed in Genomicus v64.01 (www. 
dyogen.ens.fr/genomicus-64.01). The synteny maps included all 
teleost species available in the browser. Clock genes not yet 
annotated in Genomicus were identified in Ensembl and synteny 
maps were manually constructed. 

Daily rhythm experiment 

Atlantic cod juveniles weighing 100.0 ±6.0 g (mean ± standard 
deviation [SD]) were stocked in 250 m' 5 tanks at a density of 60 
individuals per tank. To ensure minimal disturbance of the fish 
during sampling, each tank was exclusively dedicated to a single 
sampling point. Prior to the experiment, the fish were acclimated 
for at least 3 weeks under a 12L:12D daily photoperiod regime. 
Water parameters were monitored and maintained as follows: 
temperature at 7.0±0.6°C (mean ± SD) and dissolved oxygen at 
89±2.8%. A commercially formulated diet (Amber Neptun, 
Skretting AS, Stavanger, Norway) was supplied at a daily ration 
of 5% (w/w) of the fish body weight. Fluorescent white light tubes 
(Aura Light International AB, Karlskrona, Sweden) were used to 
provide illumination and light intensity was measured during a 24- 
h cycle with a Hanna Hai 97500 Luxmeter (Hanna Instruments, 
Kungsbacka, Sweden). Tissue sampling was performed at 3-hour 
intervals (Zeitgeber time: ZT0, 3, 6, 9, 12, 15, 18, 20 and 24) for a 
period of 24 h. Samples taken at ZTO were collected immediately 
after the light reached its maximum intensity (120 lux), while those 
at ZT24 were collected just before the gradual transition to the 
light phase, with an intensity increase of 2 lux-min 1 over 
60 mins. There was an approximate 30 min transition time 
(ZT12), in which the light intensity decreased at a rate of 4 
lux-min 1 to total darkness between the light (ZTO-9) and dark 
(ZT 15-24) phases. ZT12 samples were collected during this 
transition period. Ten fish were taken at each sampling point 
and killed by immersion in seawater containing 0.5 g-L 1 tricaine 
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methanesulfonate (Sigma, Oslo, Norway). Collection of samples 
during the dark phase was conducted in a room with white light 
intensity not greater than 1 lux and it was ensured that exposure of 
an anesthetized fish to these conditions did not last longer than 
5 min. Fast skeletal muscle was taken from the area below the 
second dorsal fin, skin was removed and the tissue was washed 
with cold, sterile 1 x PBS to remove contaminating blood before 
immersion in liquid nitrogen. Samples were stored at — 80°C until 
RNA extraction. Blood was also collected from the caudal vein for 
melatonin quantification. Heparinized blood samples were centri- 
fuged at 1,200 xg for 15 min at 4°C. Thereafter, the plasma was 
collected and kept at — 80°C until further analysis. 

Plasma melatonin assay 

Plasma melatonin was quantified by competitive ELISA with a 
commercially available Melatonin ELISA kit (IBL, Hamburg, 
Germany). This kit has been previously validated with fish plasma 
melatonin and had a method specificity of nearly 100% [25]. 

Tissue and early ontogeny sample collection 

Tissue samples from adult Atlantic cod previously collected in a 
sister study [23] were used to determine the expression of clock 
genes in various tissues and organs, namely blood, liver, spleen, 
stomach, mid-intestine, kidney, brain, pituitary, heart, gills, eye, 
dorsal skin, ventral skin, testis, ovary, and muscle. Fish were reared 
under natural photoperiod conditions (NL) for Bodo, Norway 
(67°N, 14°E) and sacrificed as above on week 45 (sunrise at 
approximately 8:08 AM and sunset: 15:23 PM). All tissue samples 
were collected between 9 am and 12 noon. 

Samples used for early ontogenetic expression of clock genes 
were previously collected in a separate study [26]. The fertilized 
egg defined the onset of ontogenesis and was denoted as t= 0. The 
various developmental stages were confirmed by visualization 
under an optical microscope and the following samples were 
collected: unfertilized egg, 1 -cell (5.2 hour post-fertilization at 7°C, 
hpf), 2-cell (6.8 hpf), 16-cell (13.5 hpf), oblong (31 hpf), germ ring 
(58 hpf), 50% epiboly (98 hpf), 10 somite (119 hpf), golden eye 
(225 hpf), hind gut (304 hpf), 1st feeding (346 hpf), and 20 days 
post-first feeding (828 hpf). The photoperiod regime was constant 
darkness (DD, 0L:24D) from unfertilized egg to hind gut stage and 
thereafter shifted to constant illumination (LL, 24L:0D, white light 
with an intensity of approximately 80 lux). 

Gene expression analyses 

Total RNA was extracted from the tissue samples using the 
mirVana miRNA Isolation Kit (Ambion, Oslo, Norway). The 
RNA was then quantified with a Nanodrop ND-1000 Spectro- 
photometer (Thermoscientific, CO, USA) and its quality assessed 
by denaturing electrophoresis on a 1.2% (w/v) agarose gel. cDNA 
was synthesized from 1 |0,g/ mL total RNA by QuantiTect Reverse 
Transcription kit (Qiagen, Nydalen, Sweden). 

All primers used in the expression studies were designed with 
PerlPrimer (www.perlprimer.sourceforge.net) across intron/exon 
borders to avoid amplification of contaminating genomic DNA, as 
reported [27]. Expression of the clock genes in different tissues and 
during early ontogeny was determined by semi-quantitative PCR 
(RT-PCR) using the primers given in Table SI. Thermocycling 
parameters were 94°C for 3 min, followed by 30 cycles of 30 s at 
94°C, 30 s at identified annealing temperature for each primer set 
(Table SI), and 30 s at 72°C, with a final elongation step of 72°C 
for 3 min. PCR products were visualized by electrophoresis on a 
1.0% (w/v) agarose gel and photographed with a Kodak gel 
documentation system v. 1.0.5 (Oslo, Norway). Elongation factor 1- 
alpha 1 (eefla) was used as an internal reference in tissue 



distribution analysis. On the other hand, luciferase (luc) was used 
as internal control in the expression of clock genes during early 
stages of ontogenesis because it was previously shown in zebrafish 
that most housekeeping genes are not stable during this period, 
particularly during maternal-zygotic transition [27,28]. 

Quantification of clock gene transcripts in the fast muscle 
during a daily cycle was performed by real-time PCR (qPCR) 
using the primers indicated in Table SI and the LightCycler 489 
SYBR Green I master chemistry (Roche, Basel, Switzerland) on a 
LightCycler 480 (Roche) thermocycler. Diluted cDNA was run in 
duplicate, including minus reverse transcriptase and no template 
controls. The qPCR reaction was performed according to this 
thermocycling protocol: initial denaturation at 95 °C for 15 min, 
followed by 45 cycles of 15 s at 94°C, 20 s at defined annealing 
temperature (Table SI) for each primer set and 20 s at 72°C. In 
order to calculate amplification efficiencies, a five-point standard 
curve of 2-fold dilution series were prepared from pooled cDNA. 
Cycle threshold (Cj) values were calculated with the LightCycler 
480 software using a fluorescence arbitrary value of 1 . To validate 
that the primers were amplifying the gene of interest, PCR 
products were verified by Sanger sequencing. In order to 
determine if there was an expressional correlation between clock 
and muscle-related genes (Table S2), transcript levels of selected 
muscle-related genes were also quantified in fast muscle during a 
daily cycle with the same PCR protocol. 

Three reference genes (Table S2), namely acidic ribosomal protein 
(arp), ubiquitin (ubi) and eefla were tested for their suitability for 
normalization of the experimental data using the geNorm 
algorithm, as reported [27]. Expression of target genes was 
normalized using the geometric average of arp and ubi, for they 
were the most stable on this given experimental setup. 

Statistical analyses 

Differences in clock gene expression between time points were 
analyzed with the SigmaStat statistical package (Systat software, 
London, UK). A one-way ANOVA was performed on data sets 
that passed the tests of normality and equal variance, and Tukey's 
multiple comparison test followed to delineate differences between 
time points. For data sets that did not follow a Gaussian 
distribution or did not meet the equal variance requirements, 
Kruskal-Wallis one-way ANOVA on ranks followed by Dunn's 
multiple comparison test were used instead. The level of 
significance was set at P<0.05. 

To evaluate daily rhythmicity in the expression of the clock 
genes and muscle-related genes, COSINOR analysis was per- 
formed by fitting a periodic sinusoidal function to the expression 
values of genes across the nine time points, according to the 
formula / (t) = M+ Acos(t/ pi/ 1 2 - cp), where / (t) is the gene 
expression level at a given time, mesor (M) is the mean value, A is 
the sinusoidal amplitude of oscillation, t is time in hours and (p is 
the acrophase (peak time of the approximating sinusoidal 
function). The statistical significance P of the approximated 24 h 
waveform was defined by the noise/signal of the amplitude 
(SE(A)/A) and expression was considered to display daily 
rhythmicity if P<0.3 [29]. Correlation analysis between clock 
and muscle-related gene expression was done using the Pearson's 
correlation test (r). 

Results 

Atlantic cod clock genes and their evolutionary 
characteristics 

Partial sequences of 1 8 clock genes from both loops of the clock 
machinery were cloned in Atlantic cod: amtll (bmall, Genbank 
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Figure 1. Radiation trees of representative clock genes from the transcriptional activation (A: arnth and repression (B: per) arms of 
the Atlantic cod molecular clock system. The unrooted trees were constructed by maximum likelihood analysis using an LG substitution model 
with four substitution rate categories. Branch support was calculated by aLRT SH-like tests. Orthologous genes from teleosts and tetrapods were used 
in the construction of the phylogenetic trees. Atlantic cod clock genes cloned in the present study are highlighted in red bold font. Teleostean clades 
are circled by a dotted red line. Phylogenetic trees of other gene families are provided in Figure S1A-I. 
doi:10.1371/journal.pone.0099172.g001 



accession JN643704), amtl2 (bmal2, JX035864), clock gN643707), 
cryla gN643708), crylb (KC204818), cry-dash (KC204819), cry2 
JN643709), cry3 (JX035865), npasl (KC204820), npasS 
(JX035866), nrldl (rev-erba,, KC204822), nrld2a {rev-erbjla, 
JN643712), nrld2b (rev-erbfib, KC204823), perl (JN643710), per2a 
(JN643711), per2b (JX035863), timeless (ton) JN634817) and rora 
(KC204821). 

Teleost arnth formed two distinct clades. The fish arntll group 
was closely related to their tetrapod orthologs with short branch 
lengths (Fig. 1A). Three clusters of arntl2 were observed: piscine 
arntl2a including G. morhua amtl2, piscine arntl2b and tetrapod 
arntl2. The fish arntl2a clade containing G. morhua arntl2 also 



included orthologs from Oreochromis niloticus, 0. latipes and T. 
nigroviridis. The three paralogs of Atlantic cod per genes were found 
in three clusters (Fig. IB). G. morhua perl was closely related with its 
orthologs in D. rerio, 0. latipes and O. niloticus. The two Atlantic cod 
per2 genes clustered in two separate groups, along with other 
teleost orthologs, such G. aculeatus, O. latipes, O. niloticus, T. 
nigroviridis, T. rubripes and Paralicthys olivaceus. The other clock genes 
found in Atlantic cod were also closely related to their teleost and 
other vertebrate orthologs (Fig. S1A-I). 

Comparative analysis of the genes surrounding each clock 
paralog revealed the conservation of synteny amongst teleosts, 
namely in D. rerio, G. aculeatus, 0. latipes, 0. niloticus, T. nigroviridis, T. 
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Figure 2. Partial synteny map of arnt/2 of Atlantic cod showing high syntenic conservation. Orthologous genes in G. morhua, 0. latipes, T. 
rubripes, T. nigroviridis, O. niloticus and X. maculatus, are indicated by block arrows showing their position and orientation. Additional synteny results 
can be found in Figure S2A-0. (un: unidentified gene). 
doi:10.1371/journal.pone.0099172.g002 
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Figure 3. Early ontogenetic expression of clock genes in Atlantic cod. The cDNAs used in semi-quantitative RT-PCR were from the pooled 
RNA of approximately 50 specimens from each stage [26]. Luciferase (luc) was used as an external reference. 
doi:1 0.1 371 /journal.pone.00991 72.g003 
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Figure 4. Presence of clock gene transcripts in different tissues of Atlantic cod. Expression of clock genes in various tissues (blood, liver, 
spleen, mid-intestine, kidney, brain, pituitary, heart, gills, eyes, dorsal skin, ventral skin, testis, ovary, and fast muscle) was analyzed by semi- 
quantitative RT-PCR. Eefla gene was used as an endogenous reference. 
doi:1 0.1 371 /journal.pone.00991 72.g004 
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Figure 5. Changes in the plasma melatonin level during a daily cycle. Plasma melatonin levels were quantified by competitive ELISA. Values 
shown are mean+SEM of plasma melatonin (n = 6). The background represents the photoperiod regime: the light background is the light phase, the 
dark background is the dark phase and light gray indicates the light-dark transition phase. 
doi:1 0.1 371 /journal.pone.00991 72.g005 



rubripes and Xiphophorus macuktus (Fig. S2A-I). Synteny was not 
conserved in the genomic regions containing crylb and nrld2b but 
could be clearly observed for amtl2, rora and tint. The represen- 
tative partial synteny map for amtl2 shown in Fig. 2 revealed high 
conservation, with the presence of several genes {mdt4b, nrlh4, 
scll7a8, pfibpl, si:dkey, polr3b, rfx and ric8b) in the same genomic 
location and orientation as in other teleost species. Nevertheless, 
this synteny profile was less conserved in medaka. 

Early ontogenetic expression of the clock genes 

The various molecular components of the circadian clock 
system in Atlantic cod were first detected at different ontogeny 
stages (Fig. 3). amtl2, crylb and cry2 transcripts were maternally 
transferred and detected in all stages of early development 
examined, albeit at different levels. Per2a and per2b were also 
maternal transcripts but their levels gradually decreased as 
maternal-zygotic transition progressed; these genes were detected 
again during the latter part of the development in a stage- 
dependent manner. Some clock gene transcripts were only 
detected during the latter part of the development. For instance, 
clock and rora started to be transcribed only at the golden eye stage, 
whereas arntll and cryla mRNAs were only found at 20 days post- 
first feeding. The two npas paralogs as well as cry3, nrld2a and 
nrld2b could not be detected during early ontogeny. 

Presence of clock gene transcripts in different peripheral 
tissues 

Clock gene transcripts were found in multiple tissues of adult 
Adantic cod and showed some tissue specificity in their expression, 
except amtl2, cry2 and per2a, where ubiquitous expression was 
observed in all tissues studied (Fig. 4). With the exception of npasl 
and rora, all clock genes were detected at varying levels in fast 
skeletal muscle. The majority of clock genes were highly expressed 
in brain, with the exception of cry-dash. The clock genes arntl2, 
clock, per2a, per2b and nrldl were prominent in organs related to the 
endocrine system (e.g. pituitary) and light perception (e.g. eyes). 
Apparent differences were observed in expression patterns 
between paralogs. For example, cryla was found in 13 out of the 
1 8 tissues while crylb was only present in 1 1 tissues. Similarly, per2a 
was detected in all tissues while per2b was only found in 1 1 tissues. 
Interestingly, the majority of clock genes was present at a higher 



level in dorsal rather than in ventral skin. Regardless of sex, arntll, 
clock, per2a and nrldl were detected abundandy in the gonads. On 
the other hand, cry3, perl, per2b and tint transcripts were only found 
in ovary but not in testis. 

Rhythmicity of clock gene transcription during a daily 
cycle 

Plasma melatonin levels showed a biphasic diurnal pattern, with 
higher concentrations during the dark phase than the light phase 
(Fig. 5). Its highest level was observed just after the light intensity 
reached its maximum (ZTO), while the lowest level was noted 
3 hours prior to light-dark transition (ZT9). As darkness 
progressed, the level of melatonin increased throughout the dark 
phase, reaching a significantly constant level at later hours (ZT18 - 
ZT24). 

There was no significant difference on the temporal expression 
of arntll but the arntl2 paralog displayed daily rhythmic expression 
(P— 0.01) and it was preferentially expressed during the dark phase 
having an acrophase at ZT 16.1 h. mRNA levels of clock, sl known 
co-dimer of arntl, displayed daily rhythmicity (P= 0.06) although 
with an acrophase during the light phase (ZT 9.04 h). Expression 
of npas2, an alternative component of the positive arm of the 
transcriptional feedback loop, showed daily rhythmic expression 
(P=0.08) with an acrophase at ZT 4.36 h. The rhythmicity 
parameters of these genes are provided in Table 1. Transcript 
levels of 3 out of the 9 transcriptional repressors studied (cry2 
(P=0.05), ay3 (P=0.27) and per2a (P=0.07)) displayed daily 
rhythmic expression (Fig. 6, Table 1). Cry2 had an acrophase at 
ZT 1.04 h, which occurred in an inverse phase to arntl2. 
Significant temporal differences and daily rhythmicity were 
observed in cry3 transcript levels. Per2a was preferentially expressed 
during the dark phase with the acrophase noted at ZT 17.5 h. 
Significant temporal differences were observed in crylb, tim, perl 
andper2b mRNA levels but these changes were not rhythmic. Two 
of the three orphan nuclear receptor genes exhibited daily 
rhythmic expression: nrldl (P=0.21) and nrld2a (P=0.01) 
(Fig. 6). Both genes had a light phase responsive profile with 
acrophases at ZT 6.01 h and ZT 9.01 h for nrldl and nrld2a, 
respectively. Expression oinrld2b was arrhythmic (P=0.51), even 
if some significant temporal differences were observed. 
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Figure 6. Expression of clock genes in fast skeletal muscle during a daily cycle. The values are mean+SEM (n = 6) of the normalized 
transcript levels of each clock gene. Statistical difference between time points is indicated by different letter notations. The broken line is the periodic 
sinusoidal function of gene expression in a circadian cycle constructed from the periodicity parameters calculated by COSINOR. An asterisk (*) besides 
the gene name indicates significance of daily rhythmicity. The photoperiod regime is represented by the composite block above the graph. White, 
black and gray represent the light phase, the dark phase and the light-dark transition phase, respectively. 
doi:10.1371/journal.pone.0099172.g006 



Daily expression of muscle-related genes and its 
correlation with expression of clock components 

my/5 (P=0.15) and mbnll (P= 0.05) displayed daily rhythmic 
expression (Fig. 7) with acrophases at ZT 7.5 h and ZT 7.2 h, 
respectively. Even if other muscle-related genes did not display 
daily rhythmicity, their expression profiles had expression peaks 
during the light phase, with the exception of mstn (Table S3). 

Expression of the clock components showed intra-group and 
inter-group positive/negative correlations (Table 2). Particularly, 
there was a strong positive correlation within components of the 
transcriptional repressor arm specifically in the following gene 
pairs: cry3:per2b and per2a:crylb. The components of the transcrip- 
tional activator arm arntll, amtl2 and npas2 showed moderate to 
strong negative correlation to several members (i.e. cry2, ay3, tim 
and per2d) of the transcriptional repressor arm. 

The patterns of mRNA levels for most muscle-related genes 
correlated with the daily expression of clock genes (Table 3; Table 
S4). Transcript levels of myJ5, mjhc, mbnl and foxk2 positively 
correlated with the expression of clock with myf5 and mbnll, 
showing a strong correlation index (r>0.80; Table 3). In addition, 
both daily rhythmic nuclear receptors positively correlated with 
myf5 and mbnll expression. Members of the transcriptional 
repressor arm (i.e. per2a, per2b, crylb, cry3 and tim) showed 



moderate to strong negative correlation with pena, my/6, mstn, 
myoG, mbnll and foxk2. 

Discussion 

This is the first report of the molecular components of the 
circadian clock in Atlantic cod. We have identified 18 clock genes, 
namely 5 members of the transcriptional activation arm [arntll, 
arntl2, clock, npasl and npas2), 9 genes of the transcriptional 
repression arm [cry la, crylb, cry-dash, cry2, cry3, perl, per2a, per2b and 
tim) and 4 orphan nuclear receptors of the stabilizing loop [nrldl, 
nrld2a, nrld2b and rora). These genes were transcribed in fast 
skeletal muscle, with the exception of npasl and rora, suggesting the 
existence of a possible peripheral clock in this tissue (Fig. 8). 

The circadian clock system in fish is more complex than its 
mammalian counterpart because of the existence of multiple 
copies of clock genes [1], consistent with our findings in Atlantic 
cod. The existence of 2 paralogs (i.e. cry and per genes) per 
mammalian ortholog could be explained by the teleost-specific 
third round of the whole genome duplication [30]. In contrast to 
mammals, which have only 2 copies oi cry gene [31], teleostean cry 
paralogs constitute a highly divergent family. Gene duplication is a 
major adaptive genomic response [32] and this phenomenon may 
lead to the establishment of lineage-specific traits and to the 
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Table 1. Rhythmicity parameters of clock gene transcription 
in Atlantic cod fast skeletal muscle. 





Gene 1 


Peak of expression/Acrophase (h) 


P 


arntll 


14.6 


0.41 


arntl2 


16.1 


0.01 


clock 


9.04 


0.06 


npas2 


4.36 


0.08 


cry la 


15.6 


0.62 


cry lb 


19.3 


0.63 


cry-dash 


10.3 


0.99 


cry2 


1.04 


0.05 


cry3 


2.53 


0.27 


perl 


21.2 


0.88 


per2a 


17.5 


0.07 


per2b 


21.4 


0.83 


tim 


6.33 


0.37 


nrldl 


6.01 


0.21 


nr1d2a 


9.01 


0.01 


nr1d2b 


19.3 


0.51 



'Clock genes whose transcript levels displayed daily rhythmicity are highlighted 
in bold. The P value is defined as the noise/signal ratio of the oscillation 
amplitude. P<0.3 and a significant difference with ANOVA (p<0.05) indicate 
daily rhythmicity. 

doi:1 0.1 371 /journal.pone.00991 72.t001 



development of novel biological functions [33,34]. We can 
speculate that clock genes in Atlantic cod may be sharing similar 
functions as shown by the strong evolutionary relationship of these 
clock genes with other teleost species and high syntenic 
conservation. On the other hand, differences in daily expression 
of several molecular clock components (discussed below) indicate 
some degree of divergence and functional diversification. Also, 
some paralogs (e.g. cry3, cryla, nrld2a) seem to be evolving much 



faster than their mammalian orthologs, based on our phylogenetic 
analysis results. Whether there is a redundancy or uniqueness in 
the function of each copy of the clock genes in Atlantic cod is yet to 
be ascertained in further studies. Within teleosts, the circadian 
clock system has been functionally characterized mostly in 
zebrafish [1]. The possibility that functional clocks differ across 
fish species is very likely, since the environment plays a crucial role 
in circadian regulation. 

The detection of clock gene transcripts at stages prior to 
activation of zygotic transcription indicated that some components 
of the circadian system of Atlantic cod, specifically amtl2, ay lb, 
cry2, per2a, and per2b are of maternal origin. In addition to 
maternal contribution, environmental stimuli such as light could 
influence the development of the clock system. The establishment 
of circadian rhythms clearly requires exposure to environmental 
zeitgebers such as light and temperature changes [1] and also 
depends on the developmental stage of the fish [9,35,36]. The 
present study determined the presence of clock components at a 
specific ontogeny stage but not the onset of rhythmicity of these 
genes in Adantic cod. It is noteworthy that crylb was ubiquitously 
expressed during early Atlantic cod ontogenesis. Cryptochromes 
are pterine/flavine-containing proteins and play an important role 
in phototransduction and circadian photosensitivity in DrosophUa 
[37]. crylb was prominently expressed at the golden eye stage and 
the development of this photosensitive organ may have a 
significant influence in the expression of this gene. In addition to 
the notable expression of crylb, clock and rora started to be 
expressed from the golden eye stage while perl, per2a and tim 
started to be detected after this stage only. The development of this 
photoreceptive organ could have a direct influence on the 
development of several molecular components of the circadian 
system. Other than light, food could provide a zeitgeber in the 
entrainment of circadian rhythm and this has already been 
documented in zebrafish and cavefish [14,38]. The possible 
presence of multiple peripheral clocks in Adantic cod was 
supported by the ubiquitous expression of clock genes in different 
tissues and organs. The majority of these genes were expressed 
distinctively in the eyes and to some extent in the pituitary gland, 
which was expected since these organs are related to the photo- 



Table 2. List of clock gene pairs showing either positive or negative expression correlation in fast muscle during a daily cycle. 



Positive correlation 






Negative correlation 






Gene 


Gene 


r 


Gene 


Gene 


r 


arntll 


arntl2 


0.51 


arntll 


cry2 


-0.60 


arntl2 


per2a 


0.64 


arnt/2 


npas2 


-0.78 


clock 


nrldl 


0.66 


arnt/2 


cry2 


-0.54 


clock 


nr1d2a 


0.60 


arnt/2 


cry3 


-0.73 


npas2 


cry2 


0.58 


arnt/2 


tim 


-0.60 


npas2 


cry3 


0.83 


npas2 


per2a 


-0.54 


npas2 


tim 


0.61 


npas2 


nr1d2b 


-0.57 


cry2 


cry3 


0.61 


crylb 


cry3 


-0.58 


cry3 


per2b 


0.76 


per2a 


cry3 


-0.52 


per2a 


crylb 


0.74 


cry3 


nr1d2b 


-0.74 


per2a 


nr1d2b 


0.73 









NOTE: Only correlations with r>+0.5 or r< — 0.5 and including at least one gene with significant daily rhythmicity are shown in this table. Underlined genes names 
indicate that they are part of the stability loop of the transcriptional-translational feedback loop and bold face that they are part of the positive arm. Names of genes 
involved in the negative arm are neither in bold nor underlined. The following values were set to define the degree of correlation: data are moderately correlated if 0.5< 
r<0.69 and there is a strong correlation when r >0.70. 
doi:1 0.1 371 /journal.pone.00991 72.t002 
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Table 3. Correlation of clock and muscle-related transcript levels. 





Muscle-related gene 


Clock gene 


















arntl2 


clock 


npas2 


cry2 


cry3 


per2a 


nrldla 


nr1d2a 


myf5 


-0.10 


0.84 


0.45 


0.17 


0.11 


-0.47 


0.77 


0.64 


mbnll 


-0.28 


0.82 


0.48 


-0.05 


0.05 


-0.54 


0.52 


0.55 



NOTE: Only genes displaying significant daily rhythmicity are given in this table. The overall correlation profile is given in Table S4. The arbitrary categories that define 
the degree of correlation are as follows: data are moderately correlated if 0.5<r<0.69 and the correlation is considered strong when r >0.70. 
doi:1 0.1 371 /journal.pone.00991 72.t003 




mstn foxk2 14 mbnll * a 




ZT (h) 

Figure 7. Daily expression of muscle-related genes. The values presented here are mean+SEM (n = 6) of normalized expression. Statistical 
differences between time points are indicated by different letter notations. The broken line is the periodic sinusoidal function of the gene expression 
in a daily cycle constructed from the rhythmicity parameters provided in Table S3. An asterisk (*) beside the gene name indicates that the expression 
is daily rhythmic. The photoperiod regime is represented by the composite block above the graph. White, black and gray represent the light phase, 
the dark phase and the light-dark transition phase, respectively. 
doi:1 0.1 371 /journal.pone.00991 72.g007 
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Figure 8. Molecular components of the clock system identified in fast skeletal muscle of Atlantic cod and myogenic genes with 
daily rhythmic expression. The peripheral clock components in Atlantic cod fast skeletal muscle that have been identified in this study are shown 
inside the green circle, which represents a muscle fiber. They comprise members of the transcriptional activator arm (in red: arntll, arntl2, clock and 
npas2), transcriptional repressor arm (in blue: cry?o, crylb, cry-dash, cry2, cry3, perl, per2a, per2b and tim) and the stabilizing loop (in yellow: nr1d1, 
nr1d2a and nr1d2b). Clock genes with a colored background displayed a daily rhythmic expression in fast skeletal muscle. The daily rhythmicity of 
myf5, a gene for myogenic lineage specification, and mbnll, a gene for terminal muscle differentiation, suggests a possible circadian clock control of 
myogenesis in Atlantic cod. The grey box indicates that the mechanism underlying this regulatory process remains to be identified. 
doi:10.1371/journal.pone.0099172.g008 



endocrine axis of the circadian system. The circadian clock system 
has been regarded as a fundamental mechanism in the photo- 
neuroendocrine regulation of reproduction in fish [39]. The high 
mRNA levels of clock genes such as arntll, amtl2, clock, nrldl and 
per2a observed in testes and ovaries imply that these molecular 
clock genes may influence gonadal physiology in Atlantic cod. 
Some of the clock genes identified in this study were also 
differentially regulated in gonads during a maturation cycle in 
Atlantic cod, thus supporting their possible involvement in sexual 
maturation (Lazado et al, unpublished data). The several clock 
genes identified may have immunological implications inAtlantic 
cod because they were highly expressed in immune-related organs 
and tissues, such as blood, liver, mid-intestine, gills, kidney, spleen 
and skin. In fact, it has been shown that the immune response in 
fish is regulated in a daily or seasonal pattern, and is affected by 
photoperiod conditions [40,41,42]. 

Recent studies support the notion of a decentralized clock in the 
understanding of biological rhythms. In particular, the peripheral 
clocks present in various tissues have been of great interest, as 
these oscillators are regarded as direct circadian regulators of 
genes involved in the physiology of specific tissues [17]. Circadian 
rhythms in tissues and organs are essential to ensure that 
physiological processes, such as metabolism and cell cycle, 
undergo at the optimal time as required by the organism. Clock 
genes are expressed in fast skeletal muscle in several species and 
the daily rhythmicity observed in a number of these genes provides 
a possible implication on their essential regulatory functions 
[6,19,21]. Nevertheless, the current knowledge on the relevance of 
circadian rhythm in fish skeletal muscle physiology is very limited. 
The two Atlantic cod arntl genes had a night-biased expression, 
since their peaks of expression were observed during the dark 
phase. This photosensitivity profile is corroborated by the 
expression of their orthologs in zebrafish skeletal muscle [21] but 



not in either central or peripheral clocks in rainbow trout [12] and 
mouse [43]. In the present study, there was no clear relationship 
between arntl and clock/npas2, as shown by the weak correlation 
between these transcripts and the apparent differences in their 
acrophases. A comparable expression profile was also observed in 
several peripheral tissues of zebrafish [44] . This antiphase profile 
revealed that the transcriptional activation arm of the possible 
clock system in the fast muscle is probably not similar to the 
established central pacemakers known in other fish species and 
Arntl may have another co-dimer in the activation arm of the 
clock system in fast skeletal muscle. In mammals, Bmall (Arntll) 
has the ability to interact with other bHLH-PAS factors besides 
Clock [45,46]. A similar mechanism may be present in Atlantic 
cod fast skeletal muscle but our study was limited to the 
transcriptional level. 

Among the 5 cry genes identified in Atlantic cod, their expression 
peaks indicated that cry2 and cry'i are light-biased whereas crylb is 
dark-biased. The light-biased expression of cry2 was similar in 
peripheral clocks of other teleosts, such as European seabass 
(Dkentrarchus labrax) [47]. Only cry2 and cry3 had daily rhythmic 
expression in Atlantic cod fast skeletal muscle. Both genes displayed 
moderate to strong negative correlation with arntl genes and their 
expression was in inverse phase with arntl as well. These patterns 
suggest that cry2 and cry'3 may have a significant function in this arm 
of the transcriptional feedback loop in skeletal muscle. There was 
also a moderate positive correlation between cry2 and cry3, 
suggesting that there may be a tight co-regulatory control between 
these 2 paralogs in the transcriptional repression arm. 

Cry-dash and tim are two newly-characterized clock-gene 
homologs and their involvement in the piscine clock system is 
yet to be identified. Nevertheless, cry-dash has been shown to be 
light inducible in zebrafish [48] . The dark-biased expression of the 
3 period genes was in agreement with the expression of other teleost 



PLOS ONE | www.plosone.org 



10 



June 2014 | Volume 9 | Issue 6 | e99172 



Daily Rhythm in Atlantic Cod Fast Muscle 



orthologs in central and peripheral tissues such as in goldfish 
(Carassius auratus) [29], European sea bass [10] and zebrafish 
[19,38]. On the other hand, this circadian preference is not similar 
to what has been observed in central and peripheral clocks of 
Senegalese sole [11]. In the established circadian feedback loop 
model, the proteins translated by per and cry genes interact to 
inhibit the transcriptional activation arm [1]. Per2a, the only period 
gene with daily rhythmic expression in the present study, did not 
show a positive correlation with any of the cry genes with daily 
rhythmic expression. However, the transcript levels of Atlantic cod 
per2a showed a strong positive correlation with the arrhythmically 
expressed crylb. Also, there was an inverse correlation between 
transcript levels of Atlantic cod per2b and clock/ npas2, two genes 
with circadian rhythmicity in the activation arm. In mammals, 
Clock has a constitutive expression while Bmal displays daily 
rhythmic variations in the suprachiasmatic nucleus [49]. Such a 
mechanism suggests that some genes with arrhythmic expression 
could still be important components of the peripheral clock in 
Adantic cod muscle. 

Nuclear receptors such as nrldl, nrld2 and rora are modulators 
in the stabilizing loop of the circadian mechanism [1]. Unlike a 
previous report in zebrafish [21], rora was not detected in Atlantic 
cod fast skeletal muscle but two orphan nuclear receptor genes, 
nrldl and nrld2a, had daily rhythmic expression with preference 
towards the light phase. This expression profile was inversely 
related to those of arntl paralogs, in agreement to the known 
inverse relationship between these genes [43] . clock and npas2 were 
also suggested in mammals to be regulated by these nuclear 
receptors in a repressive manner [50,51], In the present study, the 
expression profile of clock and npas2 had a similar trend with the 
nuclear receptors, and clock even showed a positive correlation 
with nrldl and nrld2a. We did not explore to what extent these 
nuclear receptors render feedback stability on the circadian clock 
transcriptional-translational network. Nevertheless, the rhythmic- 
ity observed in their expression is a possible indication of their 
essential role in the circadian rhythm of fast skeletal muscle in 
Adantic cod. 

Most of the muscle-related genes studied displayed temporal 
changes in transcript levels during the daily cycle, but only myf5 
and mbnll demonstrated daily rhythmic expression. Such changes 
imply that muscle physiology is influenced by the light-dark cycle 
and could be under the regulation of clock genes. It has been 
demonstrated that the Clock/Bmall circadian regulation of MyoD 
is critical in maintaining skeletal muscle phenotype and function in 
mice [19]. However, myoD in Atlantic cod skeletal muscle did not 
show daily rhythmicity, as observed in zebrafish [21]. It is 
plausible that the circadian regulation of piscine myoD is not under 
the same mechanism reported in mouse. Other myogenic 
regulatory factors might be playing the circadian clock-related 
transcriptional control of myogenesis in fish, such as the daily 
rhythmicity of myf6 transcript levels in zebrafish [2 1] and myf5 in 
Adantic cod. Higher transcript levels of most myogenic genes 
analyzed in this study were observed during the light phase. This 
may explain, at least partly, the known positive effect of light in 
somatic growth in Atlantic cod [22,23]. The moderate to strong 
correlations found between the expression of muscle-related and 
clock genes could be linked to a possible circadian clock-regulated 
expression of muscle genes, which has been suggested in zebrafish 
skeletal muscle [21]. In mouse, approximately 7% of the muscle 
transcriptome displays circadian rhythmicity [6] and it is very 
likely a number of muscle-related genes under circadian control 
remain to be identified in fish. 



Conclusions 

This is the first study discussing the possible relevance of a 
functional clock system in fast skeletal muscle of a teleost other 
than zebrafish. We reported that several clock genes are 
transcribed in fast skeletal muscle of Atlantic cod, suggesting the 
presence of a complex peripheral clock in this tissue (Fig. 8). It is 
plausible that muscle development and growth may be under clock 
control, as evidenced by the daily rhythmicity of my/5 and mbnll 
and strong correlation observed between transcript levels of 
muscle-related and clock genes. Taken together, our data provide 
a novel comparative perspective to our current knowledge about 
the implication of circadian clocks in muscle physiology and raises 
important questions. In particular, the autonomy of this peripheral 
clock and the extent of its importance in regulating the muscle 
transcriptome remain to be ascertained. 
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the stabilizing loop (G: nrld H: rora). The unrooted trees were 
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